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ABSTRACT 

We estimate the power spectrum of SZ(Sunyaev-Zel'dovich)-effect-induced 
temperature fluctuations on sub-degree scales by using the cross correlation be- 
tween the three-year WMAP maps and 2MASS galaxy distribution. We produced 
the SZ effect maps by hydrodynamic simulation samples of the ACDM model, 
and show that the SZ effect temperature fluctuations are highly non-Gaussian. 
The PDF of the temperature fluctuations has a long tail. More than 70% power 
of the SZ effect temperature fluctuations attributes to top ~ 1% wavelet modes 
(long tail events). On the other hand, the CMB temperature fluctuations basi- 
cally are Gaussian. Although the mean power of CMB temperature fluctuations 
on sub-degree scales is much higher than that of SZ effect map, the SZ effect 
temperature fluctuations associated with top 2MASS clusters is comparable to 
the power of CMB temperature fluctuations on the same scales. Thus, from 
noisy WMAP maps, one can have a proper estimation of the SZ effect power at 
the positions of the top 2MASS clusters. The power spectrum given by these 
top wavelet modes is useful to constrain the parameter of density fluctuations 
amplitude as- We find that the power spectrum of these top wavelet modes of 
SZ effect on sub-degree scales basically is consistent with the simulation maps 
produced with ag = 0.84. The simulation samples of cr$ = 0.74 show, however, 
significant deviation from detected SZ power spectrum. It can be ruled out with 
confidence level 99% if all other cosmological parameters are the same as that 
given by the three-year WMAP results. 

Subject headings: cosmology: theory - large-scale structure of the universe 



1. Introduction 

The thermal Sunyaev-Zel'dovich (SZ) effect is due to the inverse Compton scattering of 
the cosmic microwave background (CMB) photons by hot electrons. This effect is propor- 
tional to the pressure of the electron gases, p = ksnT, where n and T being the number 
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density and temperature of electrons. Therefore, SZ effect can be used as a measurement of 
the hot baryon gases of the universe (e.g. Birkinshaw 1999; Carlstrom et al. 2002). Since 
clusters are the hosts of high pressure gases, the angular power spectrum of the SZ effect is 
very sensitive to the amplitude of mass density fluctuations parameter a 8 (Refregier et al. 
2000; Bond et al. 2002). The amplitude of the SZ effect angular power spectrum is found 
to be proportional to crl(Qbh) 2 and is almost independent of other cosmological parameters 
(e.g. Seljak et al. 2001; Komatsu & Seljak 2002). Therefore, the amplitude of the SZ effect 
power spectrum can be used as an effective constraint on the parameter cr 8 . The SZ effect 
power spectrum on scales of a few arcminutes has been used for this issue (Mason et al.2003; 
Readhead et al. 2004; Dawson et al. 2006). 

Recently released WMAP three-year data (WMAP III) refines most results of cosmo- 
logical parameters given by the WMAP first year data. They found <r 8 = 0.74±0.05 (Spergel 
et al. 2006), which is significantly lower than a 8 = 0.84 ± 0.04 of the first year data. The 
new number of os is a challenge to the cosmological parameter determined with samples of 
galaxies and galaxy clusters, most of which yield as — 0.9 — 1 if the matter content of the 
universe fl m < 0.3 (Reiprich & Bohringer 2002; Hoekstra et al. 2002; Refregier, et. al. 2002; 
Van Waerbeke et. al. 2002; Bacon et al. 2003; Bahcall & Bode 2003; Seljak, et al. 2005; 
Viel & Haehnelt 2006; McCarthy et al. 2006; Cai et al. 2006). This problem motivates us 
to study the constraint on a 8 given by the power spectrum of the SZ effect. 

We try to extract the information of the SZ effect power spectrum on sub-degree scales 
from the WMAP III data themselves. At the first glance, it seems to be hopeless to ex- 
tract such information from noisy WMAP III data, because the power of SZ effect induced 
temperature fluctuations on sub-degree scales is much less than that of CMB temperature 
fluctuations (e.g. Cooray et al. 2004). A direct analysis on the WMAP data also shows that 
the SZ effect contribution to the CMB fluctuations on scale of the first acoustic peak should 
be less than 2% (Huffenberger et al. 2004; Spergel et al. 2006). 

This problem can be solved if considering that the SZ effect map is highly non-Gaussian, 
while CMB is Gaussian. Although the mean power of the thermal SZ effect on sub-degree 
scales is no more than 2% of the power of CMB, the SZ temperature fluctuations at local 
modes associated with clusters are comparable to the CMB temperature fluctuations. Thus, 
the SZ temperature fluctuations power on sub-degree scales would be detectable from noisy 
WMAP maps if one can identify such local modes. To implement this detection, we developed 
a DWT (discrete wavelet transform) algorithm for the power spectrum estimation. The DWT 
analysis is powerful to pick up weak but non-Gaussian signals from Gaussian background 
(Donoho 1995). 

We will use the cross- correlation between the WMAP III data and 2MASS galaxies. 
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Since the redshift of 2MASS XSC galaxies is up to z ~ 0.1, at which clusters are on angular 
scales of sub-degree. Moreover, the thermal SZ effect on angular scales of clusters has been 
significantly identified with the 2nd and 4th order cross-correlations between the first year 
WMAP data and 2 MASS XSC galaxies (Afshordi et al. 2004, 2005; Myers et al. 2004; Cao et 
al. 2006). The WMAP-2MASS cross-correlation would be effective to extract the SZ power 
on sub-degree scales from the WMAP III maps. 

The paper is organized as follow. In §2, we analyze the non-Gaussian features of the 
SZ effect maps generated with cosmological hydrodynamic simulation. §3 describes the 
identification of the SZ effect of 2MASS clusters from the WMAP-2MASS cross correlation 
with the wavelet modes. §4 first describes the algorithm of estimating the SZ power spectrum 
from the WMAP III maps, and then presents the result of the SZ effect constraint on a 8 . 
The discussion and conclusion are given in §5. 

2. The power spectrum of SZ effect fluctuations 

2.1. Maps of SZ effect 

To develop the method of estimating power spectrum, we first study the statistical 
properties of the SZ effect temperature fluctuations. We produce the SZ effect maps with 
the WIGEON (Weno for Intergalactic medium and Galaxy Evolution and formatiON) code, 
which is a cosmological hydrodynamic iV-body code based on the Weighted Essentially Non- 
Oscillatory (WENO) algorithm (Harten et al. 1986; Liu et al. 1994; Jiang & Shu 1996; Shu 
1998; Fedkiw, Sapiro & Shu 2003; Shu 2003). The details of this code can be found at Feng 
et al. (2004). Some WIGEON samples have already been used in cosmological studies (e.g. 
He et al. 2004, 2006; Zhang et al. 2006; Liu et al. 2006). 

The samples are simulated in a 100 h^ 1 Mpc cubic box with 512 3 meshes for gas, and 
equal number of particles for dark matter. We use the standard ACDM model, which is 
specified by the matter density parameter Q m = 0.27, baryon density parameter fib = 0.044, 
cosmological constant fi A = 0.73, Hubble constant h = 0.71. In order to test the effect of 
the amplitude of density fluctuations, we produce two samples with <j 8 = 0.74 and 0.84, 
respectively. The ratio of specific heats is 7 = 5/3. The transfer function is calculated 
using CMBFAST (Seljak & Zaldarriaga 1996). In order to correct the underestimate of rare 
massive galaxy clusters due to the finite box size, the initial conditions are generated by 
sampling a convolved power spectrum following the method of Pen (1997) and Sirko (2006), 
which accurately model the real-space statistics, such as a 8 . 

The atomic processes including ionization, cooling and heating are modeled as the 
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method in Theuns et al. (1998). We take a primordial composition of H and He (X=0.76, 
Y=0.2A) and use an ionizing background model of Haardt & Madau (1996). Star forma- 
tion and their feedback due to supernova explosions and AGN activities were not taken into 
account. 

Since the redshift of 2MASS galaxies is up to 0.1, we produce the SZ maps by integrating 
the pressure of gas from z — to 0.1. The simulation data is stored for every redshifts 
corresponding to comoving distance nx lOO/i^Mpc, n being integral. Three simulation boxes 
are stacked together and integrated along different angular direction to compose a SZ map. 
Two 2-D maps of the Compton parameter y for og = 0.74 and 0.84 are shown in Figure 
1. Both maps have a size of 19.1° x 19.1°, with resolution of 512x512. We produce 40 
independent maps. The total area, 40 x (19. 1°) 2 , is comparable with the area of 2MASS 
data we used (§3.1). 

Figure 2 gives the angular power spectrum for two samples with a% = 0.74 and 0.84, 
respectively. We see that the amplitude of the power spectrum is significantly dependent 
on a 8 . Comparing Figure 2 with other simulations (e.g. Refregier et al. 2000; Seljak et al. 
2001; Komatsu & Seljak 2002), our result shows about the same power as others on scales 
equal to and larger than 0°.5, while it is much lower than others on scales <C 0°.5. This 
is expected, as the power on scales less than 0°.5 mainly comes from clusters with redshift 
> 0.1. If considering all the contribution from z = 6, our simulation samples show the same 
behavior as others, i.e., declining around 0°.l. 

2.2. The DWT power spectrum of SZ effect map 

As shown by Figure 1, the value of y is very small in most space, while high value 
of y appears only in small and localized areas. One can then take the advantage of DWT 
analysis to describe the field y(x) (x = (xi,x 2 )) with DWT decomposition. The DWT 
variables of a map y(xi,x 2 ) in the area of < X\ < Li and < x 2 < L 2 are given by the 
decompositions of scaling functions 0j,i(x) and wavelet functions ^^(x). We use Daubechies 
4 wavelet (Daubechies (1992), Fang & Thews (1998). The functions 0j,i(x) and VluM can 
also be found in Yang et al. (2001). 

Briefly, The DWT scaling function, 0j i(x), is a window function for cells (modes) in 
the area L^/2^ < x x < L^k + l)/2* and L 2 l 2 /2 h < x 2 < L 2 (l 2 + l)/2- ?2 . The scale index 
j — (ji ? j2), ji and j 2 can be any positive integer. The position index 1 = (h,l 2 ), h and 
l 2 can be 0, 1... 2 jl — 1, and 0, 1... 2 J2 — 1, respectively. On the other hand, the wavelet 
function is admissible, satisfying J ^(x^x = 0. 
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For the field y(x), the DWT variables are 

yj ' 1= /^x)dx / y(x) ^ l(x)dX ' (1) 

ej, = | y(x)^,i(x)dx, (2) 

Since 0j,i(x) is a window function of cell (j, 1), y- h \ is the mean of y(x) in cell (j, 1). Because 
■0j,i( x ) is admissible, €- x measures the fluctuation of field y(x) with respect to the mean, y- J: \, 
in the mode (j,l). For a given j, the mean of (ejj) over all modes is zero (Fang & Thews 
1998). Since the basis ip^i are orthogonal and complete, the DWT variables of the field y(x) 
do not lost information and cause false correlations. 

With the DWT variables, the power spectrum of the 2-D field y(x) can be measured as 
(Pando & Fang 1998; Fang & Feng 2000): 

2-*l-l;2*2-l 

-fj = (l e j,i| 2 ) = 2ii+j2 I^ij2;ii,/ 2 | 2 ' ( 3 ) 

ii=0;i 2 =0 

|^ 1 j 2 -; 1 ; 2 | 2 is the power at mode (j,l) = (ji, j 2 ; h,h)- Therefore, the DWT power spectrum 
P } actually is the mean of local powers | ej , i | 2 over all 2 J1+J2 modes on scale j. It has been 
shown that the DWT power spectrum can be expressed by the Fourier power spectrum as 
(Fang & Feng 2000; He et al. 2005a): 

1 oo oo 

P j = ^ E E mnJ^M^Pin^n,), (4) 

n\ =— oo ri2=— oo 

where ip(n) is the Fourier transform of the basic wavelet ip(x). Because we have the nor- 
malized relation ^2 n n | "0 ( J^i / 2 J_1 ) -0 (ri2 / 2 J " 2 ) | 2 / 2- 7 " 1 +J ' 2 = 1, the DWT power spectrum P- } is 
equal to the Fourier power spectrum banded on scales around (LiL 2 /2- ?1+J2 ) 1 / 2 . Figure 3 
shows both the angular power spectrum C\ and DWT power spectrum P- } for the simulation 
samples. The angular scale of j is 19.1/2- 7 degree, and therefore, / ~ 9.4 x 2 j . Figure 3 shows 
that DWT power spectrum is the same as the angular power spectrum C t . We will use the 
DWT analysis on angular scales of sub-degree, on which the effect of spherical surface is not 
significant. 

2.3. Non-Gaussianity of SZ effect map 

The DWT power spectrum can effectively detect the non-Gaussian feature of a field 
(Jamkhedkar et al. 2003). Taking this advantage, one can reveal the non-Gaussianity of the 
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SZ effect field y(x) with the DWT power spectrum of the SZ field y(x) defined as 

top 

where summation of "top" goes over only modes having high values l^^iuhl 2 among the 
2h+h modes for a given j. In Figure 4, we plot the Pj and Pj top of the simulation samples, 
in which top power spectrum, Pj top , is given by top 1%, 0.3%, 0.1% among the 2 jl+32 modes. 
We see that all the power spectra Pj top are close to Pj. For samples with <r 8 = 0.84, there is 
more than 70%, 50%, 35% power attributes from top 1%, 0.3%, 0.1% modes, respectively. 
That is, the power |ej° p | 2 of each mode among the top 1% modes should be as high as about 
70 times of the mean (|ej i| 2 ) = Pj. These results show that the field y(x) is highly non- 
Gaussian. From Figure 4, one can see that for sample with as = 0.74, all the power spectra 
Pj top of top 1%, 0.3%, 0.1% are less close to Pj. This is because the clustering of a 8 = 0.84 
sample is stronger than that of a 8 = 0.74. 

Because (e^i) = 0, the value of (le^il 2 ) 1 / 2 actually is the standard deviation(a) of the 
probability distribution function (PDF) of ey. Thus, the events of 70(|ej 5 i| 2 ) actually are ~ 
8<T events. All the top 1% wavelet modes are 8a events. All the top 0.3% and 0.1% modes 
correspond, respectively, to ~ 170(|ej 5l | 2 ), and ~ 350(|ej 5 i| 2 ), or 13 a and 19 a events. For a 
Gaussian field, the 10 a events should be no more than 10~ 23 . Therefore, the PDF of e^i has 
a very long tail. In other words, a significant part of the SZ effect DWT power spectrum is 
given by the long tail events. These long tail events make the SZ effect on sub-degree scales 
to be detectable. Although the mean of SZ power is less than about 2% of the power of 
CMB, for long tail events, the SZ power is comparable with that of CMB. 

One can estimate the power spectrum of SZ effect temperature fluctuations with the top 
modes (long tail event) of eji. For observed maps, however, not all the top modes of are 
given by SZ effect. To reduce this uncertainty, we add a condition that the useful modes (j, 1) 
for top ej ; i should also be the modes significantly showing SZ effect by the WMAP-2MASS 
cross correlation. It is equal to assume that the power of SZ effect temperature fluctuations 
mainly comes from clusters (Komatsu & Kitayama 1999; Komatsu & Seljak 2002). We 
check this assumption with the relation between variables y- Jt \ and ej i. We found that the 
top modes of \e-^\\ are largely coincident with the top modes of the Compton parameter y- } \. 
For sample of a 8 = 0.84, 84% top modes of y- 3 ^ on the angular scale of 0°.5 are coincident 
with the top modes of \e^\\ on the same scale. This result indicates that modes with high 
i generally are also the modes of high local power of SZ effect temperature fluctuations. 
High i generally is related to clusters. 

The fact that the coincidence between top y^i and top ej,i is not 100% is due partially 
to the so-called warm-hot intergalactic medium (WHIM) (Dave et al 2001; He et al. 2004, 
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2005b). It shows that a significant fraction of hot baryon gases (T > 10 5 ) are not located in 
collapsed structures with high density contrast p/p > 100, but in the areas with a median 
overdensity 10 < p/p < 30. The WHIM is also the source of SZ effect. Nevertheless, the 
major part of top modes of \e- h \\ is associated with clusters. Therefore, a proper method to 
estimate the SZ power spectrum would be to identify the modes with top numbers of both 
y- h i and | ej,x | - The details of this method will be given in the following two sections. 

3. SZ effect identified from WMAP-2MASS correlation 

3.1. Data 

The data of galaxies in the 2MASS extended source catalog (XSC, Jarrett et al. 2000) 
used for cross-correlation analysis is the same as Guo et al. (2004) and Cao et al. (2006). 
That is, we use galaxies of K_m_k20fe, which measures the magnitude inside a elliptical 
isophote with surface brightness of 20 mag arcsec -2 in i^-band. There are approximately 
1.6 million extended objects with K s < 14.3. To avoid the contaminant of stars, we use a 
latitude cut of \b\ > 10°. We also removed a small number of bright (K s < 9) sources by 
the parameters of the XSC confusion flag (ccJlag) and visual verification score for source 
(vc). To eliminate duplicate sources and have a uniform sample, we use use_src = 1 and 
dup_src = 0. We also use a cut of 10.0 < K s < 14.0 to ensure the sample is complete 
greater than 90%. This sample contains 987,125 galaxies with redshift up to z ~ 0.1. 

For the CMB temperature fluctuations, we use the full resolution coadded 3 year sky 
map (Hinshaw et al. 2006) for W band. The maps have had the ILC(Internal Linear 
Combination) estimate of the CMB signal to highlight the foreground emissions (Hinshaw et 
al. 2006). We also use the three annual maps to construct three difference maps between the 
1st and 2nd, 2nd and 3rd, 3rd and 1st year, which are useful to check statistical significance. 

As in Cao et al. (2006), we subject both the 2MASS and WMAP maps to an equal-area 
projection. We then have 2-D maps of galaxy number density, n(x), and CMB temperature 
fluctuations, AT(x) with coordinate x = (xi, x^). As has been shown, the equal-area projec- 
tion will not affect statistical analysis on sub-degree scales, but helpful for using the DWT 
method. We analyze the cross correlation between n(x) and AT(x) in two 123°. 88 x 123°. 88 
areas, which are, respectively, in the northern and southern sky. The 2MASS XSC galaxies 



1 The notations of the 2MASS parameters used in this paragraph are from the list shown in the 2MASS 
Web site http://www.ipac.caltech.edu/2mass/releases/allsky/doc 
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are resolved to 10 arcsec. The beam size of WMAP W band map is 0°.22g, while the scale 
of galaxy clusters at the median redshift of the 2MASS samples is 0°.5. 

In order to apply the method of DWT power spectrum (§2.2), we also use the DWT vari- 
ables to describe the maps n(x) and AT(x). They are wavelet variables e?j = J n(x)^ji(x)<ix, 
e^i = f AT(x)'0j ) i(x)dx, and scaling function variables n- h \ = f ri(x)0j i i(x)(ix/ f 0j i i(x)c?x, 
ATj 5 i = / AT(x)0j,i(x)dx/ / 0j,i(x)dx. The mode (j, 1) occupy a area of 123°.88/2 J1 x 
123 .88/2J' 2 at position around [Zi(123°.88)/2 il , Z 2 123°.88)/2^]. The angular distance be- 
tween modes 1 and 1' at scale j is given by 9 = 123.88°|1 — l'|/2- 7 . 

The variables e? y and e| j describe, respectively, the fluctuations of the CMB temperature 
and galaxy number density in the mode (j,l), while the variables ATj i and n- h i describe, 
respectively, the mean temperature and the mean number density of galaxies in the mode 
(j,l). Therefore, n- 3j \ and ATj i can be seen as the maps of galaxies and CMB temperature 
functions smoothed on scale j. 



3.2. DWT clusters 

We picked up top clusters on scale j by the top members of n- h \. These clusters are called 
DWT clusters. The details of identification of clusters have been studied with simulation 
and real samples (Xu et al. 1999, 2000). It showed that the clusters identified by DWT 
n- } \ are statistically the same as the clusters identified by the friend-of-friend method if the 
mean size of the friend-of-friend identified clusters is the same as DWT clusters. The clusters 
identified by the friend-of-friend method usually have very irregular shapes (Jing & Fang 
1994) and are inconvenient to estimate the statistical significance of the cross correlation. 
On the other hand, the DWT variables of 2MASS galaxies n- h \ and WMAP map AT^ are in 
the same cell, the statistical significance of the n^i-ATy cross correlation is unambiguous. 
Scaling functions are orthogonal with each other, different DWT clusters consist of different 
galaxies. This is also necessary for statistical analysis which we will take below. With this 
method, we identify top clusters on scales j = (8,8), corresponding to scale ~ 1.8 h _1 Mpc 
of the 2MASS samples. 



2 http: / /lambda. gsfc.nasa.gov/product 
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3.3. Thermal SZ effect of 2MASS clusters 

The SZ effect of 2MASS galaxies can be measured with the cross-correlation between 
the WMAP and 2MASS clusters defined as Cao et al.(2006): 

AT wmap-2massQ 1 _ ^ = (CtfATy , V ) , (6) 

where the variable Cy\ is given by: 

f 1 if (j,l) is a cell of identified clusters, , . 

j '' ~~ \ otherwise. ^ ' 

The average (...) covers all the DWT clusters on scale j. Therefore, AT wmap - 2mass (|l - l'|) 
measures the mean of the CMB temperature fluctuations on a distance |1 — l'| from identified 
DWT clusters of 2MASS data. When |1 - l'| > 1, AT wmap - 2mass (|l - l'|) actually is the mean 
of AT on a large area, and therefore, it should approach to zero. Figure 5 presents the cross 
correlation AT wmap - 2mass (|l - l'|) between top 100 j = (8, 8) DWT clusters and WMAP map, 
in which the scales j' of ATjy/ are taken to be j' = (8, 8), j' = (8, 7), and j' = (7, 7) from top 
to bottom. In all case ATjy is from the full resolution coadded 3 year WMAP map of W 
band. The error bars are given by the 1-cr standard deviation of the 100 CMB simulations. 

Figure 5 shows typical anti-correlation of SZ effect. At the position of DWT clusters, 
i.e. |1 — 1'| = 0, the mean of temperature change ATyi/ is always negative, which is expected 
from the thermal SZ effect on W band. For j' = (8, 8), AT wmap - 2mass (|l - l'| = 0) ~ -25±15 
/iK, while for j' = (8,7) and j' = (7,7), the temperature decreases are —20 ± 16 fiK and 
— 18 ± 14 fiK respectively. These results are consistent with that given by cross-correlation of 
the first year WMAP map and 2MASS galaxies (Myers et al. 2004; Hernandez-Monteagudo 
et al. 2004; Cao et al. 2006). 

To further check the SZ signals, we calculate the cross-correlation of eq.(6) by using 
the maps given by the difference between the 1st and 2nd, 2nd and 3rd, 3rd and 1st year. 
The results are shown in Figure 6. All the anti-correlations shown in Figures 5 completely 
disappear in Figure 6. It confirms that the temperature decreases at the positions of the 
DWT clusters in the three annual maps are about the same and not due to noise. 

To view the clusters richness dependence, Figure 7 plots AT wmap_2mass (|l — l'| = 0) 
against the number of top clusters, n. As expected, the higher the richness, the stronger the 
SZ effect. AT wmap ~ 2mass (|l — l'| = 0) rapidly decreases with the number of top clusters, and 
approaches to the noise level at n > 300. We should, however, point out that the positive 
correlation between the clusters richness and SZ effect is statistical. That is, the order of 
clusters richness is not, one-by-one, corresponding to the order of the amount of SZ effect 
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temperature decrease. Consequently, the relation of richness-SZ effect in the range n < 100 
is not very stable. When n > 100, the relation of richness-SZ effect, as shown in Figure 6 
is very stable. Therefore, in all the statistics below, the number of top clusters should be 
> 100. 

In the DWT analysis, wavelet modes are orthogonal from each other. Therefore, the 
area fraction of top 300 j = (8, 8) DWT clusters is 300/(2 x 256 2 ) = 0.23%. As shown in §2, 
the mean power of SZ effect of top 0.3% modes would be 170 times higher than mean power 
of total mode. On the other hand, on sub-degree scales, the CMB power is higher than the 
SZ effect power by a factor ~ 10 2 . Therefore, we can expect that the mean of SZ temperature 
decreases of the top 0.3% modes should be larger than the CMB temperature fluctuations 
by a factor ~ a/170/100 ~ 1.3. Figure 8 plots the richness-dependence of the S/N ratio, i.e. 
the ratio between SZ effect temperature decrease and the 1 — a error bar shown in Figure 7. 
Since the error bar of Figure 7 is given by the simulation of CMB temperature fluctuations, 
and therefore, it actually is to measure the the CMB temperature fluctuations. We see from 
figure 8 that the S/N ratio is about 1.3 - 1.6. This result is consistent with the estimation 
mentioned above. 



4. Estimation of DWT power spectrum of SZ effect 

4.1. Method 

The observed CMB temperature fluctuations are given by: 

AT(x) = AT cmb (x) + AT sz (x) + AT sccond (x), (8) 

where AT cmb (n) is the primeval temperature fluctuations, AT sz (x) is the SZ effect and 
AT sccond (n) is due to secondary effects other than the SZ effect, such as the ISW effect and 
radio point sources. Since the random fields AT cmb (x), AT sz (x) and AT sccond (x) basically 
are uncorrelated from each other, the DWT power spectrum, Pj, of the observed AT(x) 
consists of three terms as: 

p. = p. cmb + p. sz + p. second , (9) 

where 

2M-1;2*2-1 

(10) 



p>cmb,sz, second ^ \ r ~cmb,sz,second 

i 1= o-i 2 =o 



Since p. cmb is about 10 2 times larger than P? z on angular scale of ~ 0°.5, one cannot 
extract P? z from noisy map AT(x) if the field AT sz (x) is Gaussian. As shown in §2, the SZ 
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effect DWT power spectrum is largely given by the top modes. The local power of the top 
DWT modes can be as large as 10 2 times of their average. Therefore, the SZ effect local 
power of the top modes is comparable to i^ cmb . On the other hand, the primeval temperature 
fluctuations AT cmb (x) are Gaussian, and the PDF of local temperature fluctuations ej\ is 
Gaussian with zero mean and standard deviation (/y 1110 ) 1 / 2 ^ rp^g p r i meV al temperature 
fluctuations map does not have the modes with very high local power. Thus, the local power 
of modes associated with top clusters would be measurable from the noisy map of AT(x). 

The measurability of power of top modes can be further tested as follows. Subjecting 
eq.(8) to a wavelet transform, we have: 

~cmb | ~sz | -second / 1 i \ 

e j,l = e j,l + e j,l + e j,l > ( U ) 

For most modes (j, 1) on scales of clusters, \e?J b \ is much larger than |e^j|, and therefore, 
the latter is negligible. But for modes corresponding to the top clusters, |ey| would be 
comparable to |ej™ b |. On these modes, the term |e?*j cond | is also negligible, as statistically the 
positions of secondary effect are not coincident to top clusters. Thus, on these modes, we 
have: 

I ~ 1 2 I ~cmb 1 2 , I ~sz 1 2 /i<)\ 

l e j,l| - l e j,l I + l e j,ll » 

and therefore, 

^— Y id' 2 = jr- Y i^i 2 - w- Y l g j'i nb l 2 ' ( 13 ) 

JVt °P top JVt °P top iVt °P top 

where the summation "top" is over modes with top eji and also the top DWT clusters. iV t0 p 
is the number of top modes. Because CMB field is statistically isotropic (uniform), Gaussian 
and uncorrelated with 2MASS data, the average of the primeval temperature fluctuations 
power lej™ 13 ! 2 done by (l/N top ) J2to P should be about the same as the average over all modes. 
We have then: 

El^l^^opdej,,! 2 ), (14) 
top 

where { | e j 5 1 1 2 ) is the mean of |ej 5 i| 2 over all modes on scale j. Thus, considering le?^ 2 > 0, 
eq.(13) yields: 



Ei^i 2> <i^i 2 >' ( 15 ) 



Eq.(15) requires that the mean local power over modes at top clusters should be larger than 
the mean on all modes. 

We test Eq.(15) by the cross correlation defined similar to eq.(6): 

Tj#-11) = ^## , (16) 
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where C- Jt \ is given by eq.(7). The numerator on the r.h.s. of eq.(16) is to measure the 
mean local power of CMB temperature fluctuations over the modes of DWT clusters on 
scale j'. The denominator of eq.(16) is for normalization. The inequality in eq.(15) requires 
Tj/(|1 — l'|) > 1 at |1 — l'| =0, i.e. the local power of the modes, which shows significant SZ 
temperature decrease, should also show power larger than the mean (le^i'l 2 ) over the whole 
sky. The result of Tj/(|1 — l'|) is plotted in Figure 9, in which the scale j' of local power |ej',i'| 2 
is taken to be j' = (8, 8), j' = (8, 7) and j' = (7, 7). For each scale, we count 0.11% modes 
with top |ej',i'| 2 , and all the counted modes are located in the positions of the 300 top 2MASS 
clusters, corresponding to 150, 75, and 38 top modes of the scales j' = (8, 8), j' = (8, 7) and 
j' = (7,7). Figure 9 shows that all the mean powers of CMB temperature fluctuations on 
the three scales are significantly larger than the mean (lejyl 2 ). This power excess is due to 
the SZ effect of the 2MASS clusters, which are larger than the noise of WMAP III data. 

4.2. Results 

With these preparation, we can estimate the SZ effect power spectrum given by modes 
associated with the top DWT clusters of 2MASS galaxies. It is, 

top 

where the "top" is taken the same modes as shown in Figure 9. Using eqs.(13) and (14), 
eq.(17) gives, 

(18) 

This is the basic formula of extracting the SZ temperature fluctuations power spectrum of 
2MASS clusters from the WMAP III map. The result is presented in Figure 10, in which 
we show the Pj top given by 0.11% top modes(black points). The scales of j = (8,8), (8,7) 
and (7,7) is corresponding to angular 0°.48, 0°.68 and 0°.96, respectively. The error bars are 
given by jack-knife. As a comparison, Figure 10 also shows the simulation results of both Pj 
(dashed lines) and Pj top for top 0.11% modes(solid lines). The error bars are given by the 40 
independent simulation maps using jack-knife. 

The shape of all the power spectra shown in Figure 10 are similar. Their amplitudes 
are, however, very different. The extracted Pj top of SZ effect from the 0.11% top modes 
associated with top 2MASS DWT clusters is much higher than Pj top of the 0.11% top modes 
of the simulation map with cr 8 = 0.74. The data points are about the same, and even higher 
than Pj of sample <r 8 = 0.74. Since Pj is given by all modes of the sample, while Pj top is only 



ptop 



1 



top 



ml 
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from the top 0.11% mode. From eqs.(3) and (17) we must have Pj > Pj° P . Therefore, the 
result of sample a 8 = 0.74 is unacceptable with confidence level 99%. On the other hand, for 
sample of a 8 = 0.84, the extracted SZ power spectrum (data points) is safely lower than the 
simulation Pj, and it is also basically consistent with 0.11% top modes Pj top of the simulation 
map. 

Therefore, we can conclude that the power spectrum of the SZ temperature fluctuations 
on sub-degree scales extracted from the WMAP III shows a significant deviation from sample 
with the parameter of density fluctuations amplitude <t 8 = 0.74, but is consistent with the 
simulation sample produced with erg = 0.84, if all other cosmological parameters are the 
same as that given by the three-year WMAP. 

5. Discussions and conclusions 

We studied the SZ-effect-induced temperature decrease and fluctuations of the three- 
year WMAP maps caused by the 2MASS galaxies. The DWT variables are very useful to 
study this problem, as the two sets of DWT variables, given by the decompositions of the 
scaling function and wavelet, are just suitable to describe, respectively, the SZ temperature 
decrease and temperature fluctuations. 

We show that the maps of both SZ temperature decrease and temperature fluctuations 
are highly non-Gaussian. The PDF of the DWT local power has a long tail. More than 
70% power of the SZ effect on clusters scales is given by 1% top modes associated with top 
clusters. On the other hand, the field of CMB temperature fluctuations is Gaussian. One 
can measure the long tail events of the SZ effect from the noisy WMAP III maps. These 
long tail modes yield an useful estimation of SZ effect power spectrum. It gives, at least, a 
robust lower limit of the power spectrum of the SZ temperature fluctuations. 

We found that the WMAP temperature shows significant decrease at the positions from 
100 up to ~ 300 top clusters of 2MASS galaxies. We also found that there is significant 
power excess at the positions of these top clusters. This power excess is a measurement of 
the SZ effect power spectrum. We find that the power spectrum given by the 0.11% top 
modes associated with top clusters of the 2MASS galaxies, shows even higher amplitude 
than that of simulation sample with cr 8 = 0.74. On the other hand, the power spectrum of 
these top modes is bascially consistent with cr 8 = 0.84 simulation SZ effect map. Therefore, 
the SZ effect power spectrum on sub-degree scales seems to favor parameter <r 8 > 0.74. This 
result supports the direct measurement of the SZ effect power excess on arcminute scales 
from data of CBI (Mason et al.2003; Readhead et al. 2004) and BIMA (Dawson et al. 2006). 
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Fig. 1. — The SZ effect Compton parameter y maps of simulation samples with redshift 
z < 0.1 and a 8 = 0.74 (top) and a 8 = 0.84(bottom). 
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Fig. 2. — Angular power spectra of SZ effect y maps for samples given by hydro dynamic 
simulations with parameter a 8 = 0.74(solid) and 0.84(dashed). 
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Fig. 3. — A comparison between the angular power spectra of SZ effect parameter y maps 
given by the Fourier modes (top panel) and DWT modes (bottom panel). In each panel, the 
solid and dashed lines represent o~s = 0.74 and 0.84 respectively. I at top panel is the index 
of spherical harmonic function. Scale j at the bottom panel corresponds to angular scale of 

19.172*. 
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Fig. 4. — DWT power spectrum of SZ effect for simulation samples of <7g = 0.74(top) and 
0.84 (bottom). In each panel, we show the DWT power spectra of whole modes (solid), 
and 1% (dot-dashed), 0.3%(dashed) and 0.1% (dotted) modes with top power. Scale j is 
corresponding to angular scale of 19.1°/2 J '. 
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Fig. 5.— Cross correlation AT wmap - 2mass (|l-l'|)[eq.(6)] between the top 100 j = (8,8) DWT 
clusters of 2MASS galaxies and ATj/ i given by the full resolution coadded 3 year sky map 
for W band. The sacle of AT y>l is taken to be j' = (8,8) (top), (8,7) (middle) and (7,7) 
(bottom). In each panel, the error bars are given by the 1-a error of cross correlation between 
the top 2MASS DWT clusters with 100 CMB simulations. The angular scale of |1 - l'| is 
|1- l'|123.88°/2 8 degree. 
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Fig. 6.— Cross correlation AT wmap - 2mass (|l-l'|)[eq.(6)] between the top 100 j = (8,8) DWT 
clusters of 2MASS galaxies and ATyy given by map of the differences between 1st and 2nd 
years (solid), 2nd and 3rd years (dotted), 3rd and 1st years (dashed) of WMAP data on W 
band. The sacle of ATj/ i is taken to be j' = (8, 8) (top), (8, 7) (middle) and (7, 7) (bottom). 
The angular scale of |1 - l'| is |1 - l'|123.88°/2 8 . 
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Fig. 7. — The cluster richness dependence of the SZ effect temperature decrease 
AT wmap - 2mass (0), in which n is the number of top j = (8,8) clusters of 2MASS galaxies. 
The scale of ATy^ is taken to be j' = (8,8) (top), (8,7) (middle) and (7,7) (bottom). The 
error bars are given by the 1-a error of cross correlation between the top 2MASS DWT 
clusters with 100 CMB simulations. 
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Fig. 8. — The signal-to-noise ratio of the SZ effect temperature decrease AT wmap_2mass (0) 
against clusters richness, in which n is the number of top j = (8, 8) clusters of 2MASS 
galaxies. The scale of ATy^ is taken to be j' = (8,8) (top), (8,7) (middle) and (7,7) 
(bottom). 
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Fig. 9. — Tj/(|1 — l'|) [eq.(16)]. The average is taken over 0.11% modes with top |ej',i'| 2 in 
the positions of 300 top cluster with j = (8, 8). The scales j' is taken to be j' = (8, 8) (top), 
j' = (8, 7) (middle) and j' = (7, 7) (bottom) . The angular scale of |1 - l'| is |1 - l'|123.88 /2 8 . 
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Fig. 10. — The DWT power spectrum of SZ effect temperature fluctuations at W band. 
The results from eq.(18) with 0.11% top modes associated with the top DWT clusters is 
shown by black point with error bar. The dashed and solid lines show the power spectrum 
of simulation sample with all modes and 0.11% modes respectively. The top panel is for 
sample of cig = 0.74, the bottom is for <7g = 0.84. Scale j is corresponding to angular scale of 
123.88 /2 J . The three black points of scales of j = (8,8), (8,7) and (7,7) have angular scale 
of 0°.48, 0°.68 and 0°.96, respectively. 



